A GEOMETRIC VIEW OF THE RATIONAL LANDEN 
TRANSFORMATIONS 

JOHN HUBBARD AND VICTOR MOLL 



Abstract. We provide a geometric interpretation of a new rational Landen 
transformation and establish the convergence of its iterates. 



1. Introduction 

The transformation theory of elliptic integrals was initiated by Landen in [6, 7], 
wherein he proved the invariance of the function 

GM = / I 2 2a ,2 • 2. M 

Jo vo 2 cos 2 # + 6 2 sm 6 



under the transformation 



Gauss [4] rediscovered this invariance in the process of calculating the arclength of 
a lemniscate. The limit of the sequence (a n , b n ) defined by iteration of (1.2) is the 
celebrated arithmetic-geometric mean AGM(a, b) of a and b. The invariance of the 
elliptic integral (1.1) leads to 

= G(a,b). (1.3) 



2 AGM(a, b) 



General information about AGM and its applications is given in [3]. A geometric 
interpretation of the transformation (1.2) is given in [5]. 

A transformation analogous to the Gauss-Landen map (1.2) has been given in 
[1] for the rational integral 

U 6 {a 1 ,a 2 ; b ,b 1 ,b 2 ) = / -g- ^ TVT dz - ( 1A > 
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Indeed, the integral Ue is invariant under the transformation 

(1) a\ai + 5ai + . 
h ~ (oi + a 2 + : 

(2) _ ai+a 2 + 6 



(i) Qia 2 + 5ai + 5a 2 + 9 

01 ~ (a 1+ a 2 + 2)4/3 



& 2 



2 (ai + a 2 + 2) 2 /3 

(i) 6 + 6i + fr 2 

~ (ai + a 2 + 2) 2 / 3 
b (a 2 + 2) - 

Oi 

(l) _ b + b 2 

(ai + a 2 + 2)2/3 ■ 



(i) _ b (a 2 + 2) + 2bi + 6 2 (ai+3) 
1 ~ oi + a 2 + 2 



This transformation was obtained by a sequence of elementary changes of vari- 
ables and the convergence of 

(a n ,b„) := (a[ n \ai n \b n \b{ n \b^). 

For any initial data (ao, bo) € x M.^_ there exists a number L, depending upon 
the initial condition, such that 

(a„,b„) — (3,3,L,2L,L), (1.6) 

so that 

f/ 6 (a„,b„) — > Lx \- ( 1 - 7 ) 

The invariancc of [/q under (1.5) shows that 

C/ 6 (a ,b ) = Lx|. (1.8) 

Therefore the iteration given above becomes an iterative form of evaluating the 
integral. 

The main result of [2], quoted below, is an extension of (1.5) for an even inte- 
grand. 

Theorem 1. Let R{z) = P(z)/Q(z) with 

p-i p 

P( z ) = ^6^2(p-l-j) and Q( z ) = ^ a . z 2(p-j)_ ( L9 ) 
j=0 j=0 

Define dj = for j > p, bj = for j > p — 1, 

j 

fc=0 

for0<k<p-l, 

1 p 

dl = 2^>P- fe ' (LU) 

2p-l 



= a 3 b P-l-3 + k (1-12) 



k=0 
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for < j < 2p — 1, and 



a p {i) 



Let 



.l + ELi d fe *f i = 0- 



at = 



2 2j Q(1) 2 ( 1 -Vp) 



for 1 < i < p — 1, and 

b+ = Q(lf l 'P +l 'P- 2 x 



-p-l-i , 

(cfe +c 2p -i-fc)( 

. fc=0 ^ 



p — 1 — k + i 
2i 



for < i < p — 1 . Finally, define the polynomials 



p-i 



P+(z) - Yl K^ 1 ^ and Q+{z) = ]T a+ z 2 ^ . 



k=0 



Then 



o Q(«) 



P+(£) 

o Q+CO 



k=0 



dz. 



(1.13) 



(1.14) 



(1.15) 



(1.16) 



(1.17) 



The proofs in [1, 2] are elementary but lack a proper geometric interpretation. 
In particular, the proof of (1.6) given in [1] could not be extended even for degree 
8 in view of the formidable algebraic difficulties. The goal of this paper is to show 
that the transformation (1.14, 1.15) is a particular case of a general construction: 
the direct image of a meromorphic 1-form under a rational map. This will allow us 
to prove an analogue of (1.6,1.8) for the integral 



U2p{ 



a,b) := f 
Jo 



00 b Q z 2 P- 2 + hz 2 P~ 



z 2 p + aiz 2 P- 2 H hi 

In fact, we prove that the sequence x„ starting at 

x = (ai, • • • , a p -i; bo, ■ ■ ■ , b p -\) 
and defined by x„ + i = x+ satisfies 



■ dz. 



(1.18) 



P 
p-1 



where 



- U 2p (a, b). 

7T 



p-1 
p-1 



Moreover the convergence of the iteration is equivalent to the convergence of the 
initial integral. 
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2. The direct image of a 1-form 

Let 7r : X — > Y be a proper analytic mapping of Ricmann surfaces (i.e., a finite 
ramified covering space), and <p be a tensor of any type on X. Then 7r*<£ is the 
tensor of the same type on Y, defined as follows: Let U C Y be a simply connected 
subset of Y containing no critical value of n, and let <ti, • • • , a k '■ U — > X be the 
distinct sections of 7r. Then the direct image of ir^p is defined by 



i=l 

This defines 7r»<p except at the ramification values of w, where n*(p may acquire 
poles even if <p is holomorphic. 

We will be applying this construction in the case where ip is a holomorphic 
1-form, and in this case ir^ip is analytic. 

Lemma 1. If ir : X — > Y is proper and analytic as above, and ip is an analytic 
1-form on X , then ir^ip is an analytic 1-form on Y . Furthermore, for any oriented 
rectifiable curve 7 on Y , we have 

TT*<fi= ip. 

Proof. The only problem is to show that Tr*ip is holomorphic at the critical values. 
It is clearly enough to show that the contribution of a neighborhood of a single 
critical point is holomorphic. Thus we may assume that n(z) — w = z m for some 
m, and that 

p = (a k z k + a k+ iz k+1 + . . . )dz, 

with k > 0. 

For i = 0, . . . ,m — 1 set o~i{w) = CVo(w), where C = e 2m / m and cro(w) = w 1 /™ 
for some branch of the 1/m power, for instance the one where the argument is 
between and 2n/m. Then 

. u , . I if k + 1 is not divisible by m . . 

nJz dz) = < ,, , , s , 2.2) 
v ; \u^ k+1 - m ^ m du if fe + 1 is divisible by m. X ' 

Thus the first term of the power series for ip to contribute anything to ix^p is the 
term of degree m — 1 , and it contributes to the constant term; similarly, the terms 
of degree 2m — 1, 3m — 1, . . . contribute to the terms of degree 1,2,..., all positive 
powers. □ 

This has a useful corollary. Recall that the degree of a meromorphic function 
is the maximum of the degrees of the numerator and the denominator when the 
rational function is written in reduced form. 

Lemma 2. If it : P 1 — > P 1 is analytic, and ip = R{z)dz is a meromorphic 1-form on 
P 1 so that R is a rational function of degree k, then ir^ip can be written as R\{z)dz, 
where R\ is a rational function of degree at most k. 

Proof. By Lemma 1, the number of poles of ir*ip is at most equal to the number of 
poles of ip, and clearly the orders of the poles cannot increase either. □ 

Note. It is quite possible for the degree of -n^p to be less than the degree of 
ip. This can happen in two ways: we might have poles at two points 21,22 such 
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that 7r(zi) = 77(^2), and then the polar parts at these points could cancel. We may 
also have a pole of order > 1 at a critical point, and then the order of the pole 
at the corresponding critical value will be less (in fact, the pole might disappear 
altogether) . 

3. A PARTICULAR BRANCHED COVER 

We will be concerned with the specific map 

tt(z) = ™: =f ^- (3-1) 

This mapping can also be viewed as the Newton map associated to the equation 
z 2 + 1 = 0. As such it has ±i as superattractive fixed points, and tt is conjugate to 
F{z) = z 2 via the Mobius transformation M(z) = (z + i)/(z — i); indeed Moio 
M- 1 = F. 

Let us list some properties of n. 

Lemma 3. If (p has no poles onlcP 1 , then 

/OO />OC 
ip = / n*(p. 
-oo J — oo 

Proof. If ip has no poles on M (including at infinity), then the integral converges. 
Since n maps the real axis (including oo) to itself as a double cover, the result 
follows from Lemma 1. □ 

Let r : P 1 — > P 1 be the map z \— ► —z. Then clearly 7r o r = r o n. Call ip even if 
r*ip = ip, and odd if r*<p = —<p. 

Note. When ip = R(z) dz with R a rational function, then ip is even if and only if 
R is odd, and <p is odd if and only if R is even, since dz is odd. 

Lemma 4. We have the following identities: 
(a) 

ir*ir*p = ip + r*(p. 

(b) If ip is even, then ir^ip = 0. 

(c) If ip is odd, then ir^ip is also odd. 

Thus we can restrict our attention to odd 1-forms. Below we calculate 7r* (R(z) dz), 
where R{z) is an even rational function. We will only consider the case when the 
numerator of R has degree at least 2 less than the denominator, as this avoids a 
pole at infinity which would prevent the integral over R from converging. 

The explicit evaluations of the form ir^ip described below were conducted using 
Mathematica. The corresponding sections are 



a±(w) = w± \/w 2 + 1, (3.2) 
so that for ip = &(z)dz we have 

Tr.p = <P(a + (w))^ + <P(a_(w))^-. (3.3) 
The calculations are formidable. 
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Example 1. Let 



Then 



<p = b " dz. (3.4) 
a z z + oi 



n*ip = T — 7 i \2 dw - 3 - 5 ) 

Acioaiw 2 + (a + axY 

Observe that the new 1-form can be written as 

A(a ,ai) 

7T*<£ = b X r 2 2 , r dw, (3.6) 

G^(a ,ai)w 2 + A 2 (a ,ai) 

where A(o, b) and G(a,b) are the arithmetic and geometric means of a and 6 re- 
spectively. 



Example 2. The form 



is transformed into 



b z 2 + h 

ip = -. - dz (3.7) 

a z 4 + aiz z + a 2 



8(a 2 6 + a 6i)w 2 + 2(a + ai + a 2 )(6o + &i) ^ 
16a a 2 u; 4 + 4(a ai + 4aoa 2 + aia 2 )w 2 + (ao + ai + a 2 ) 2 



4. The convergence of (71-*)™^ 

In this section we present the principal theorem of the paper. 

Theorem 2. Let ip be a 1-form, holomorphic on a neighborhood U of M C . 
Then 

dz 



lim (7r*)> = - ( / <p ) r 



w/iere i/ie convergence is uniform on compact subsets of U. 

Proof. We find it convenient to prove this for the map F(z) — z 2 , which is conjugate 
to 7r. In that form, the statement to be proved is that if ip is analytic in some 
neighborhood U of the unit circle, then 

Iim(F.)>=-Lf/ 

n^oo 27TI \Jsi / Z 

Any such 1-form </? can be developed in a Laurent series 

/ oo 

E 

where X)fcli(l a fcl + l a -fc|)p fe < oo for some p > 1. Note that 

2tti J s i 

In this form it is very easy to compute F*tp. 



k \ dz 

= | 2^ "- : 1 — • 

\k— — oo / 
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Lemma 5. The mapping on 1- forms is given by 

°° dz 

F * ( P = ^ a 2k z k — . 
^ z 

k— — oo 

Proof. This is what was computed in Equation 2.2. □ 

Thus in the "basis" of forms z k ^, the vector corresponding to k = is an 
eigenvector with eigenvalue 1, and the rest of the space is nilpotent: 

{F*) m z k — = 

z 

if m is greater than the greatest power of 2 which divides k. This comes close to 
proving Theorem 2, but it doesn't quite; for instance 

(v^ fe | dz dz 
h z )~ = 4^) 

is also fixed under F* . We cannot argue merely in terms of formal Laurent series: 
convergence must be taken into account. 

But this is not too hard. Consider the region Ur defined by < \z\ < R, and 
the space Ar of analytic 1-forms 

<p= ( J2 a k z k \ ^ 

\k— — oo / 

on IIr such that 

oo 

\\<f>\\ = \a \+Y,(\ a k\+a-k\)R k < oo. 

Wc then have 



k=l 



dz 

\K<P-ao — \ 

z 



: 5^(|a 2 »fc| + |a- 2 " fc |)i? fc 
k=i 

OO jjfc 

= E(i^ fe i + i«- 2 " fc i)i? 2 " fc ^ 

k=l 

< Jrlkll. 

This certainly shows that 7r"</?— a^^f tends to 0, in fact very fast: it superconverges 
to 0. 

□ 



5. Normalization of the integrands 

In the previous section we have produced a map 7r* of 1-forms <f> = R(z) dz that 
does not increase the degree and the integral over [0, oo]. Moreover, we have seen 
that the integrands n^ip converge as n tends to infinity. This doesn't quite tell us 
about the convergence of the coefficients of R, because of possible common factors 
and cancellations. Here we normalize the rational functions so that 7T* induces a 
convergent iteration on the coefficients. 
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We will write the integrands so that their denominators are monic and with con- 
stant term equal to 1. The latter can be achieved by factoring the constant term 
out while the former is obtained by a change of variable of the form z^Az, with 
an appropriate A. 

Example 1. For rational functions of degree 2 we obtain 



dz = I gMgo+gi) dw . {5A) 



/o a z 2 +ai Jo 4a ciiw 2 + (a + ai) 2 

This is an identity: both sides normalize to 



VoooT Jo x + 1 
Example 2. The quartic case yields 



60 *f ^ ( ^ 



b z 2 + b 1 f°° b^w 2 + frW 

dz = / — m - m dw, (5.3) 



/ ao^ + ai^ + a 2 J aj^u; 4 + ; w 2 + a 

where 



C = 8(a 2 6 + ao&i) (5.4) 

b ( p = 2(a + a 1 +a 2 )(bo + b 1 ) 

= 16a a 2 

= 4(a ai + 4a a 2 + aia 2 ) 

a 2 ^ = (a + ai+a 2 ) 2 . 



The normalization shows that the integral 

coo i, 1/2 2 i J, I/ 2 

r b a 2 ' z z + bia Q ' ^ 
Jo z 4 + a~ Q 1 ^ 2 aia 2 1 ^ 2 z 2 + 1 

equals 

(ao + ai +a 2 )~ 1/2 x 
{a 2 bo + a h)w 2 + (b + b x )aj' 1 a 2 ' 1 



L 



o w 4 + [a ai + 4a a 2 + aia 2 )a 1/,2 a 2 1 ^ 2 (a + ai + a 2 ) 1 ]w 2 + 1 

s identity can be verified dir 

f°° dx _ f°° _ 

J x 4 + 2aa; 2 + 1 ~ J x 4 + 2ax 2 + 1 ~ 2 3 / 2 v / ^TT 



dw. 



Naturally this identity can be verified directly using 

dx f°° x 2 dx 7r 



Example 3. In the case of degree 6 we obtain 
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/ 

JO 



o aoz 



DC 



6 + aiz 4 + a 2 z 2 + as 



6 z 4 + 6iz 2 + 6 2 



dz = 



f 

Jo 



r 4 1 V + 6< 1 V + 4 1) 

o a^w 6 + o^V + a^w 2 + a { 3 



(i) 



dw 



(5.5) 



where 




32(a 3 6 + a b 2 ) 

8(a 2 6o + 3a 3 6 + «o^i + «3^i + 3a & 2 + aifo) 

2(a + ai + a 2 + a 3 )(6 + &i + 6 2 ) 

64o o 3 

16(a a 2 + 6a a 3 + aia 3 ) 

4(a ai + 4a a 2 + aio 2 + 9a a 3 + 4aia 3 + a 2 a 3 ) 
(a + ai + a 2 + a 3 ) 2 . 



(5.6) 



The normalization of (5.5) yields (1.5). 
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